A Common Origin of Neutralino Stars and 
Supermassive Black Holes 



V. I. Dokuchaev* Yu. N. EroshenW 

Institute for Nuclear Research of the Russian Academy of Sciences, Moscow 



Abstract 

To account for the microlensing events observed in the Galactic halo, 
Gurevich, Zybin, and Sirota have proposed a model of gravitationally 
bound, noncompact objects with masses of ~ (0.01 -j- 1)M@. These 
objects are formed in the expanding Universe from adiabatic density 
perturbations and consist of weakly interacting particles of dark matter, 
for example, neutralinos. They assumed the perturbation spectrum on 
some small scale to have a distinct peak. We show that the existence 
of this peak would inevitably give rise to a large number of primordial 
black holes (PBHs) with masses of ~ 1O 5 M0 at the radiation-dominated 
evolutionary stage of the Universe. Constraints on the coefficient of 
nonlinear contraction and on the compactness parameter of noncompact 
objects were derived from constraints on the PBH number density. We 
show that noncompact objects can serve as gravitational lenses only at 
a large PBH formation threshold, <5 C > 0.5, or if noncompact objects 
are formed from entropic density perturbations. 



1 INTRODUCTION 



Dark objects with masses of ~ (0.01 -v- 1)M Q were detected in the Galactic 
halo when the microlensing of stars from the Large Magellanic Cloud was ob- 
served. Brown and white dwarfs and Jupiter-like planets were proposed as 
possible baryonic candidates. According to the latest data [1], dark objects 
must account for about 20% of the hidden mass in the Galactic halo. How- 
ever, the theory of primordial nucleosynthesis imposes much more stringent 
constraints on the number of such baryonic objects. Gravitationally bound, 
noncompact objects, which are also called neutralino stars, were considered 
among nonbaryonic candidates; these objects can explain some microlensing 
events with evidence of gravitational lenses being noncompact [2, 3]. The 
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hypothetical noncompact objects are the lightest objects in the hierarchical 
large-scale structure of the Universe. In the model of Gurevich, Zybin, and 
Sirota [3], these are formed immediately after the passage of the Universe to 
the dust stage. 

To reconcile the parameters of noncompact objects with data for the ob- 
served microlensing events, Gurevich et al. [3] assumed the cosmological den- 
sity perturbation spectrum on some small scale to increase sharply with an 
rms fluctuation of the order of 1 at the peak. For adiabatic perturbations, 
rms fluctuations of the order of 0.05 correspond to these fluctuations at the 
radiation-dominated stage on the cosmological horizon scale. As was shown in 
[4, 5], fluctuations with such an rms value give rise to primordial black holes 
(PBHs) at the radiation-dominated evolutionary stage of the Universe. As we 
will see below, the mass of the forming PBHs exceeds the mass of noncom- 
pact objects by several orders of magnitude. The large difference between the 
masses of noncompact objects and PBHs stems from the fact that the mass of 
the radiation contained in a fixed comoving volume at the radiation-dominated 
stage is much larger than the mass of cold dark matter (CDM) in the same 
volume at the dust stage. 

If the power spectrum of primordial cosmological perturbations is a power 
law with an index n > 1, then PBHs are formed in a wide range of masses. If, 
however, the spectrum has a peak on some scale, then PBHs are formed mostly 
in a narrow range of masses, near the mass that corresponds to this peak. A 
sharp peak emerges in the fluctuation spectrum if the inflationary potential 
V(4>) has a flat segment [6, 7]. Indeed, if the derivative at some value of the 
scalar field is V = dV{4>)/d(j) — > 0, then a peak emerges in the perturbation 
spectrum on the corresponding scale, because the mean density perturbation 
on the horizon scale is <5h ~ Mp 3 V 3 ^ 2 /V , where Mpi is the Planck mass. A 
similar effect is achieved in inflationary models with several scalar fields [8, 
9]. In this case, the spectrum outside the peak can have an ordinary shape, 
for example, it can be a Harrison-Zel'dovich spectrum, and can give rise to 
galaxies, their clusters and superclusters according to standard scenarios. 

Certain evidence for a deviation of the initial perturbation spectrum from a 
simple power-law shape, more specifically, for a spectral break near large scales, 
k ~ (0.06 — 0.6)h Mpc -1 , was obtained in the counting of APM galaxies and in 
the Boomerang and Maxima cosmic microwave background (CMB) anisotropy 
measurements [10]. Therefore, there is reason to suggest that the spectrum 
may also exhibit features on small scales. 

Here, we show that if the perturbation spectrum has a peak, then there 
is a clear relationship between the masses of noncompact objects and PBHs; 
a PBH mass of the order of ~ 1O 5 M corresponds to a noncompact object 
mass of the order of ~ 0.1M Q . For noncompact objects to be able to serve 
as gravitational lenses, they must originate from sufficiently large dark matter 
density fluctuations. At the radiation-dominated stage, these fluctuations log- 
arithmically increase with time and become nonlinear even at this stage. To 
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study the evolution of the fluctuations at the radiation-dominated stage, we 
use the nonlinear model proposed by Kolb and Tkachev [11] and specify the 
initial data for this model according to the linear theory. 

Based on standard astrophysical constraints on the PBH number density 
in the Universe, we obtained constraints on the rms fluctuations at the peak. 
In turn, the constraints on the fluctuations give constraints on the radius of 
a noncompact object. As was shown in [3], a noncompact object can serve 
as a gravitational lens and can produce observable microlensing events only 
if its radius exceeds the Einstein radius for this object by no more than a 
factor of 10. Stringent constraints on the coefficient of nonlinear contraction 
for noncompact objects follow from this condition. 

The inferred relationship between noncompact objects and PBHs holds only 
for adiabatic cosmological density perturbations. If the density perturbations 
are entropic, then even if there are large fluctuations in the dust component, 
the radiation density on the horizon scale is almost uniform and no PBHs are 
formed. 

It should be emphasized that when talking about a common origin of non- 
compact objects and PBHs, we have in mind not the relationship between 
individual noncompact objects and PBHs but the fact that the fluctuations 
from which the entire collection of noncompact objects and the entire collec- 
tion of PBHs originate have a common spectrum. If a PBH emerged at some 
point in space, then a noncompact object can no longer emerge at this point. 
Conversely, if there is a noncompact object, then no PBH was previously (at 
the radiation-dominated stage) formed at this point. 



2 THE FORMATION OF PBHs 

The PBH formation criterion was analytically derived by Carr [5] and con- 
firmed by numerical calculations [12, 13]. The region of space with a density 
p > p c = 3H 2 /8nG can be roughly considered to be part of the closed Uni- 
verse. Gravitational collapse of this region and the formation of a PBH take 
place if the relative radiation density fluctuation 5u = (p — p c )/p c at the time 
it goes under the horizon satisfies the conditions 

S c <5 n < 1, (1) 

where 8 C = 1/3. The left-hand inequality implies that the radius of the per- 
turbed region at the time t its expansion stops exceeds the Jeans radius ct/y/3, 
while the right-hand inequality corresponds to the formation of a PBH rather 
than an isolated universe. The mass of the forming PBH in this model is 

Mbu = ||, (2) 
where Mh is the mass within the horizon. 
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In recent years, numerical experiments have revealed the so-called critical 
gravitational collapse, during which the mass of the forming PBH is [14, 15] 



M BH = AM H (6 U - S c 



(3) 



where A ~ 3, 7 ~ 0.36, and 5 C ~ (0.65 -j- 0.7). The mass (|3p can be much 
smaller than Mh. However, as shown in [16], the PBH mass distribution 
for critical gravitational collapse is concentrated near Mbh ~ Mh, and the 
contribution of low masses to the cosmological PBH density is modest. Here, 
we consider various cases where S c lies within the range 1/3 < 5 C < 0.7. 

Following [6-8], we assume that there is a sharp peak in the fluctuation 
spectrum on some fixed (in comoving coordinates) scale £ = r/a(t). Since 
the PBH formation threshold is large, 5h > 1/3, most £-scale fluctuations 
do not collapse into a PBH but are preserved until the passage to the dust 
stage if the dark matter particle free streaming length is moderately large [3]. 
The mass Mh within the horizon as a function of the mass M x of the dust 
component in fluctuations of the same comoving scale can be calculated by 
using the standard Friedmann equations. Noncompact objects are formed on 
time scales t < t eq , where t eq is the time of equality between the matter and 
radiation densities. At this epoch, a flat cosmological model serves as a good 
approximation: 
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where 77* 2 = 2irGp eq al q /3c 2 , a eq is the scale factor at time t 
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(5) 



l+z eq = 2.32-10 4 fi m /i 2 , 



Po = 1.9-10~ 29 fi m /i 2 g sm -3 is the present cosmological 
matter density, and h is the Hubble constant (in units of 100 km s ~ 1 Mpc _1 ). 
We perform our calculations for two cosmological models: the model with 
the present matter density parameter Q m = 0.3 and the cosmological term 
Q\ = 1 — Q m ~ 0.7 and the model without the A term and with Q m = 1. The 
presence of the A term reduces only to a change in p eq ; it does not affect the 
formation of noncompact objects, because the A term contributes negligibly 
to the total cosmological density on time scales t < t eq . 
For Mh and M x , we have the expressions 



M H 



Air 



4tt 



(6) 



3 -p H (a(r/ H )0 3 , M x = yPo(aM0 3 - 

On the horizon scale, a(^n)^ = 2ctn with r/u ?p and pn = 3/3271(7^. The 
present density is po = Peq(a e q/a-o) 3 - 
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Using (H), we obtain 



Mi 



H 




(7) 



GM H ( M H 
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Denote the rms density fluctuation 5n by Ah = (<5 H ) 1//2 . The fraction of 
the radiation mass that transformed into PBHs at time £h is then [5, 8] 



d6n t 5 H 
P — — exp' 



27tAh ^ v 2A H 



erf — = — erf 



(9) 



5 C 



V2A H Vv^Ah, 



Ah , SI 



where erf(x) is the error integral. Since, according to [16], the fraction of the 
collapsing mass of the Universe for critical gravitational collapse is 0.8/3, Eq. 
@ for critical collapse is also a good approximation. 

Using (0) and (P), we can calculate the cosmological PBH density parameter 
f^BR" at the present time to: 



(3 ( 3 \V6 



BH 



22/3 V2W Mi /3 GV2p^ 6 



(10) 



3.45- ID-fJ^V f^r^r^exp 



O.lM y V0.3/ V - 6 / J c 1 2A H 



To within a factor of order unity, expression (|7|) can be derived from a 
simple estimate M H ~ M x a(t e(l ) / a(t R ) ~ M E (t eq /t H ) 1//2 , where £ H = GM H /c 3 
and t eq ~ 6 - 10 10 s is the completion time of the radiation-dominated stage. 
In the same way, we can derive ([II]) from the estimate f^BH — (3Q>{t eq ) / a,(tn) . 

PBHs are formed in the tail of the Gaussian fluctuation distribution, while 
most noncompact objects are formed from rms fluctuations. Therefore, we 
repeat that not each fluctuation, by any means, from which a noncompact 
object formed could collapse into a PBH at time £h- 
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3 THE EVOLUTION OF PERTURBATIONS 



The evolution of a radiation density perturbation at the radiation-dominated 
stage follows the law [17] 



5 r = xf(x) + 



3x 2 



x 2 + 6 dx 



4-m 



fin 



where for the growing mode f(x) = A in ji(x / V3) , j\ is the spherical Bessel 
function, A in is the normalization constant, x = krj, and k is the comoving 
perturbation wave vector. The physical perturbation wavelength satisfies the 
relations 



(12) 



On the horizon scale, A p h/2 ~ 2ct, in — tt, and we obtain from (11) 8 T = <5h 



A- m (fi, where ~ 0.817. 



For adiabatic perturbations, the perturbation in nonrelativistic matter at 
x <C 1 is 5 = 3<5 r /4. In [17], an analytic solution was found for 5 in the entire 
interval from i C 1 to a; > 1 at the radiation-dominated stage. At x ^> 1, 
this solution is 
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where 7e — 1/2 ~ 0.077 and A in = <5h/0 is the same as in ([□]). 

The applicability of (|T3"D is limited to a linear stage, i<l. When passing 
to a nonlinear stage, we will use the spherical model from [11]. In this model, 
the evolution of adiabatic perturbations is described by the equation 



1 3 
1 + -y 
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(14) 



where y = a(r])/a eq and the following parametrization was introduced for the 
radius of the perturbed region: 



r = a(r])b(r])£, 



(15) 



Here, £ is the comoving coordinate and b(r]) allows for the deceleration of 
cosmological expansion in the region of enhanced density. The quantity b in 



([TBI) can be expressed in terms of S as 

b=(l + 5)- 1/3 . 



(16) 



This relation means the passage from the Eulerian description of the pertur- 
bation evolution (TT3|) to its Lagrangian description (|T4T). 

To solve (|14D requires specifying 5; at some initial y^ b\ = (1 H- 1 ^ 3 
according to (16), and the expansion rate db/dy. In [11], Eq.flUp was solved for 
entropic perturbations when the initial velocity may be disregarded, db/dy ~ 0. 
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In our case of adiabatic perturbations, the initial velocity is large; we specify 
it according to the solution fll"3|). At x ^> 1 and y <C 1, we have 

7t / 3 \ 1/6 yc 



x 



22/3 \2nJ Mi /3 GV2pV 6 



n \- 2/3 ( h \~ 4/3 ( M \~ 1/3 



(17) 



.0.3/ V - 6 / \ M © 
We obtain from flTJ), (16), and (17) 

db 
dy 



We solve Eq. (|14]) numerically. The time y\ must be chosen in such a way that 
([13]) and (|TJ]) are simultaneously valid. An optimum choice is the time when 
Si = 0.2. A test shows that the results of our calculations change by no more 
than 15% for a different choice of 5i in the range 0.1 — 0.3. Having specified 
5i, we obtain x- x and y x from (p~3|) and (|T7|). The evolution of 5 = b~ 3 — 1 is 
illustrated in Fig. 1. 

The cosmo logical expansion of noncompact objects stops when dr/dt = 0. 
According to [11], this is equivalent to 

? = — (19) 
dy y 

Denote b and y at the time the expansion stops by 6 max and y max , respectively. 
The CDM density in noncompact objects at the same time is 

Pmax — PcqZ/max^max 

(20) 

and, consequently, the radius of the noncompact object at the stoppage time 
is 

( 3M, \ 1/3 , . 

tfmax = — — • (21) 

\4vrp max / 

The numerically calculated dependence of R max on 5 H is shown in Fig. 2. 

Thus, we have shown how the radius of a noncompact object at the time its 
cosmological expansion stops can be determined for the specified perturbation 
#h on the horizon scale. 

4 THE PARAMETERS OF 
GRAVITATIONAL LENSES 

After its cosmological expansion stops, a noncompact object begins to contract, 
with its final radius being 
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where k is the coefficient of nonlinear contraction. It is generally believed 
[18] that after the cosmological expansion of a noncompact object stops, it 
is virialized, radially contracting by a factor of 2, i.e., k = 0.5. According 
to the theory of gravitational instability, k ~ 0.3 in a multiflow region [3]. 
At present, the physical processes that could cause a noncompact object to 
contract to k < 0.3 are unknown. Therefore, we take k ~ 0.3 as the lower 
limit. 

For a noncompact object to be able to serve as a gravitational lens produc- 
ing observable microlensing events in the Galactic halo, its radius should not 
significantly exceed the Einstein radius 



R E = 2^/GM x d/c 2 , (23) 
where d ~ 20 kpc for microlensing in the halo. We define e as 

e = 10f (24) 

According to calculations [3], the inequality e > 1 must hold. If e < 1, then 
the theory comes into conflict with observational data on the light curves [3]. 
One of the microlensing events with a lens mass M x ~ 0.02M Q was studied 
in detail in [19]. A comparison of the observed and calculated light curves 
showed that if a noncompact object has no baryonic core at its center, then it 
must have the size R x = 1.6 • 10 13 cm and a compactness parameter (in our 
notation) e ~ 19. In the presence of a baryonic core with a mass of 0.05Me, 
R x = 5.7 • 10 13 cm and e ~ 4.8. It was noted in [19] that the model of a point- 
like lens for this event is also acceptable because of the large observational 
errors. 



We obtain from (0), (H), and 



kb = 10*L = L9 - 1Ql5 cm (MA (_*_) V2 r 25 ) 

#max #max \M & ) V 20k P C / 

On the other hand, using our calculations (see the preceding section), we 
derived a relationship between -R max and <5h; the quantity Sn defines the present 
cosmological PBH density according to Eq. (|T0"D. It should be noted that 
noncompact objects are formed from rms fluctuations. Therefore, 5h from the 
preceding section should be set equal to the rms fluctuation <5h = A H . Using 
our numerical calculations, we derive the dependence of £7bh on the product 
ne (see Figs. 3 and 4). 

There are several astrophysical constraints on the mass and number of 
PBHs. It follows from a constraint on the age of the Universe that ^2bh < 1- 
If PBHs provide the dominant part of dark matter in the Galaxy, then they 
must tidally interact with globular clusters by disrupting them. The PBH 
mass was constrained for this case in [20], Mbh < 1O 4 M . At ~ 1, PBHs 
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are capable of distorting the CMB spectrum if they are formed about 1 s after 
the annihilation of e + e~- pairs [5]. Mass accretion by black holes at the pre- 
galactic and present epochs contributes to the background radiation in different 
wavelength ranges. However, calculations strongly depend on the model and 
yield tt BU < 1(T 3 4 KT 1 for M BH ~ 1O 5 M . In [22], the constraint fi B H < 0.1 
on intergalactic PBHs was obtained from the condition for the absence of 
reliable gamma-ray-burst lensing events for 10 5 M Q < Mbh < 1O 9 M . A more 
stringent lensing constraint, f^BH < 0.01 for the mass range 1O 6 M < Mbh < 
10 8 M Q , was obtained from VLBI observations of compact radio sources [23]. 

Let us first consider the microlensing event studied in [19] by assuming that 
e = 0.3. If the noncompact object has no baryonic core (the vertical line ke ~ 
5.7 in Figs. 3 and 4 corresponds to this case), then the constraint ^bh < 0.1 
rules out the interpretation of this event as microlensing by noncompact objects 
for all 5 C = 1/3 4- 0.7 in the two cosmological models under consideration. In 
the presence of a baryonic core (the vertical line ke ~ 1.44), the case with 
5 C = 0.7 remains admissible in the cosmological model without the A term 
(Fig. 4). 

Consider the less stringent condition ke > 0.3, which is satisfied at e > 1 
and k > 0.3. At ~ 0.7 (Fig. 3), the constraint f2 B H < 0.1 rules out all 
models with S c < 0.5. At J1 B h > 10~ 6 , the case with S c ~ 0.7 remain possible. 
If fl\ = 0, then the constraint Qbu < 0.1 leaves a narrow region (see Fig. 4) 
with M x > O.1M and 5 C = 1/3. The constraint £7bh < 10~ 3 completely rules 
out the models with S c = 1/3, but allows the cases with 5 C > 1/3. 

5 CONCLUSIONS 

We have shown that the formation of noncompact dark-matter objects (neu- 
tralino stars) proposed in [2, 3] to account for the observed microlensing events 
in the Galactic halo must be preceded by the formation of a PBH with a mass 
of the order of ~ 1O 5 M . These PBHs and neutralino stars are formed from 
the same peak in the primordial fluctuation spectrum. 

Astrophysical constraints on the number of PBHs in the Universe allowed 
us to constrain the coefficient of nonlinear contraction and compactness pa- 
rameter for a noncompact object ke; our constraints are shown in Figs. 3 and 
4. The most stringent constraints are obtained in the presence of a cosmologi- 
cal term Q\ ~ 0.7. In this case, noncompact objects can serve as gravitational 
lenses only at a large PBH formation threshold, 5 C > 0.5, which was calculated 
in the model of critical gravitational collapse. At smaller 5 C , the model of non- 
compact objects as microlensing objects is ruled out. The constraint on ke is 
significantly relaxed if the A term is small (Fig. 4). However, this possibility 
is currently considered to be unlikely. 

It is important to note that to avoid the situation where ln(f2 B H) <^ 
requires an accurate adjustment of the parameters for the inflationary model 
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that ensures that Ah is in a narrow interval, ~ 0.05 -j- 0.12. Therefore, if the 
observed microlensing events are actually produced by noncompact objects, 
then these objects most likely originate from entropic density perturbations. 
In this case, our constraints are removed, as discussed in the Introduction. 

The constraints can also be significantly relaxed if for some reason, PBHs 
are formed in smaller quantities than that given by expression (|j). This is pos- 
sible, for example, in the case of great importance of the nonlinear effects that 
accompany the generation of metric perturbations at the inflationary stage if 
these effects suppress the appearance of large fluctuations [24]. Note, how- 
ever, that the role of nonlinear effects is presently not completely understood, 
and the results of calculations depend on the specific inflationary model. For 
example, it was found in [25] that nonlinear effects cause the PBH formation 
probability to increase, which is directly opposite to the result from [24]. 

Conversely, if the noncompact nature of lenses will be proven in the future, 
then for adiabatic density perturbations, this leads us to conclude that a large 
number of PBHs with masses of the order of ~ 1O 5 M can be formed at 
the radiation-dominated sage of the Universe. These PBHs must affect the 
evolution of galaxies and their nuclei. It may well be that one of such PBHs 
was found by the Chandra space X-ray observatory in the galaxy M 82 [26]. 
In a separate paper [27], we develop a model for the formation of galaxies 
through multiple mergers of protogalaxies with condensation centers in the 
form of massive PBHs. The mergers of galaxies and the growth of central 
massive black holes at the galactic nuclei take place simultaneously with the 
formation of galaxies. The recently found correlations between the masses 
of central black holes and the bulge velocity dispersion have been shown to 
naturally arise in this scenario. 

Since noncompact objects and PBHs originate from a common perturba- 
tion spectrum, we can in principle reconstruct the shape of the perturbation 
spectrum and determine the PBH mass function from the distribution of non- 
compact objects in mass and radius using a Press-Schechter-type theory [28]. 
Unfortunately, only a few objects were detected by microlensing, and such 
a calculation will become possible only in the future when the statistics will 
improve. 

We may consider a situation that, in a sense, is reverse to the situation 
described previously. According to the hypothesis [7], the dark halo objects 
responsible for microlensing are PBHs with masses of the order of ~ (0.01 -r- 
1)M Q . If elementary particles with a free streaming length <C £ comprise the 
remaining part of the dark matter, then miniclusters will be formed from these 
particles at the dust stage. Using formula (0), we can immediately estimate 
the minicluster mass as ~ 10~ n -v- 1O~ 8 M . Such masses are possible if the 
mass of the dark- matter particles exceeds 1 GeV [3] . The formation of PBHs 
with f^BH ~ 1 requires rms fluctuations Ah — 0.06. At time t e q, fluctuations 
in the dust component 5 ~ 1, which give rise to miniclusters, correspond to 
them. 
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FIGURE CAPTIONS 



Fig. 1. The evolution of a CDM density perturbation 5. The plot corre- 
sponds to the parameters 5u = 0.04, M x = O.1M , and Q m = 0.3. The curve 
was obtained from formula (|13|) before the point (y; = 3.96 • 10~ 6 , Si = 0.2), 



and a numerical solution of Eq. ( I4"D was used at y > y- x . The cosmological 
expansion of a noncompact object stops at the radiation-dominated stage at 
y = a/a eq ~ 0.49. The dashed line indicates the evolution of 5 according to 
the linear theory (|13|). 

Fig. 2. Radius -R max of noncompact objects at the time their expansion 
stops versus perturbation magnitude <5h on the horizon scale. Solid curves 1, 
2, and 3 correspond to the masses of noncompact objects M x = 1M , O.1M , 
and O.O1M in the cosmological model with Q m = 0.3. Dashed lines 4, 5, and 
6 were constructed for f2 m = 1 at the same masses. 

Fig. 3. PBH density parameter ^bh versus nonlinear contraction coefficient 
and compactness parameter for noncompact objects ke in the model with 
Q m = 0.3. Curves (1, 2, 3) correspond to the masses of noncompact objects 
M x = (0.01, 0.1, 1)M at S c = 1/3. Curves (4, 5, 6) correspond to 5 C = 1/2 and 
curves (7, 8, 9) correspond to 5 C = 0.7 for the same masses. The horizontal 
and vertical lines correspond to the upper observational limits Qbr = 0.1, 
^bh = 10~ 3 and to ke = 0.3, 1.44, respectively. The admissible regions lie 
below the horizontal lines and to the right from the vertical lines. 

Fig. 4. Same as Fig. 3 for the cosmological model with Q m = 1. The 
vertical lines correspond to ke = 0.3, 1.44 and 5.7. 
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